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Abstract 
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invariant functional relationship among the gauge and Yukawa couplings of the third 
generation in the symmetric phase. In the case of the minimal model, we find that the 
predicted values for the top and bottom quark masses are in agreement with the present 
experimental data for a wide range of supersymmetry breaking scales. We also find that an 
experimental accuracy of less than 1% for the top quark mass could test the corresponding 
prediction of the Gauge- Yukawa unified model. 
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1 Introduction 



The remarkable success of the standard model (SM) suggests that we have at hand a 
highly non-trivial part of a more fundamental theory for elementary particle physics. 
Since, however, the SM contains many independent parameters, it has been a challenge 
to understand the plethora of these free parameters. In Grand Unified Theories (GUTs) 
fl], 1^, the gauge interactions of the SM are unified at a certain energy scale Mgut, and 
consequently its gauge couplings are related with each other. Also the Yukawa couplings 
can be related among themselves to a certain extent. These relations among the couplings 
can yield testable predictions for GUTs P, ^ . 

However, GUTs can not relate the gauge and Yukawa couplings with each other. In 
order to achieve Gauge- Yukawa Unification (GYU), within the assumption that all the 
particles appearing in a field theory model are elementary, one has to consider extended 
supersymmetry Unfortunately, it is extremely difficult to construct a realistic model 
based on the extended supersymmetry, because the model has a real structure with respect 
to SU{2)l X U{1)y HI- In superstrings and composite models such relations, in principle, 
also exist. However, in both cases there exist open difficult problems which among others 
are related to the lack of realistic models. 

Recently, an alternative way to achieve unification of couplings has been pro- 

posed; it is based on the fact that within the framework of a renormalizable field theory, 
one can find renormalization group (RG) invariant relations among parameters which 
can improve the calculability and the predictive power of the theory. This idea is called 
sometimes the principle of reduction couplings . In this paper we would like to consider 
50 (10) supersymmetric GUTs along the lines of this unification idea. We note that all 
realistic supersymmetric 5*0(10) models have to be asymptotically-nonfree, because one 
needs a certain set of Higgses to break 50(10) down to SU{?,)c x SU{2)l x U{1)y [i IH. 

The common wisdom is that the asymptotically-nonfree theories develop a Landau pole 
at a high energy scale, a fact which inevitably suggests that the theory is trivial, unless 
new physics is entering before the couplings blow up. However, there exist arguments 
leading to a different view point; the theory converges to a well-defined ultraviolet fixed 
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point, a "new phase", instead of blowing up []. Non-abelian gauge theories could have 
the same behavior, and it might be that an asymptotically nonfree, non-abelian gauge 
theory with matter couplings can change after the critical value of the couplings its phase 
due to a certain self-adjustment of the couplings and become a well defined finite theory 
Q. In this way, a dynamical unification of couplings ||T^ can be achieved, since to 
enter into the new phase the couplings are supposed to satisfy a definite relation. It 
is natural to assume that at scales below the critical value, that is, in the symmetric 
phase of a GUT these relations among the couplings are RG invariant as a remnant of 
the dynamical unification of couplings. When the GUT enters in its spontaneous broken 
phase, these RG invariant relations serve just as boundary conditions at Mqut on the 
evolution of couplings for scales below it. Since the principle of reduction of couplings is 
based on RG invariant relations among couplings, the GYU based on this principle could 
be a consequence of the dynamical unification of couplings described above. This is a 
speculation, of course, because so far there exists no reliable and decisive calculation on 
the behavior of asymptotically-nonfree, non-abelian gauge theories ^ 

There have been recently various phenomenological studies on 5*0(10) supersymmet- 
ric GUTs without GYU [^-[^, where the top quark mass was calculated from the 
requirement of the correct bottom-tau hierarchy (recall that the top Yukawa coupling 
contributes significantly to the RG evolution of the bottom and tau Yukawa couplings). 
Unfortunately, there exists a wide range of the predicted values of Mt (160 — 200 GeV 
pop . One of the main reasons is that the bottom-tau hierarchy is experimentally known 
only with a large uncertainty (~ 10% ||23|). Suppose that the bottom-tau hierarchy would 
be precisely known and the calculated top mass would exactly agree with the experimen- 
tal value. Then there should exist a unique Yukawa coupling of the third generation in 
^See ref. [|l3j, which contain also earUer references on ultraviolet fixed points. 

^Similar phenomenon has been observed in asymptotically free theories, in which the couplings have 
to be related with each other in order for the theories to be asymptotically free and hence well defined 



in the ultraviolet limit 16|. 

■^This analysis was motivated by the observation [p^ that there might exist an infrared fixed point in 

asymptotically free QCD. 

"'Even the case of QED has not been completely clarified [T5| . 
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5*0(10) GUTs, which is consistent with experimental data. It is, however, clear that this 
Yukawa coupling can not be predicted within the conventional GUT scheme. 

With a GYU the model obtains a more predictive power, and the Yukawa couplings 
become calculable, as has been experienced in our recent studies on other models p|-pT|. 
Although the GYU proposed there is a gradual, conservative extension of the usual GUT 
scheme, it has turned out to yield successful predictions. We emphasize that this success 
is not just a consequence of the infrared behavior of the Yukawa couplings ||2^. Although 
the infrared behavior is an important ingredient for the successful predictions, we should 
stress also the significance of the field content of the theories as well as their interactions 
above the unification scale. The reason is that they contribute to the j3 functions in the 
symmetric phase which fix the structure of the GYU based on the principle of reduction of 
couplings and consequently the boundary conditions for the evolution of couplings below 
the GUT scale. Therefore, it is absolutely nontrivial that (1) there exist a unique Yukawa 
coupling of the third generation in a supersymmetric 5*0(10) GUT that is consistent with 
the experimental data, and (2) this Yukawa coupling can be calculated by means of the 
reduction of couplings. 

In this paper we will consider two different models; the first one with the Higgs super- 
multiplets of 1, 10, 16, 16, 45, 54 and the three fermion generations in 16, and the second 
one with 126 + 126 instead of the singlet which provides the Georgi-Jarlskog mechanism 
p5| , pr2| to obtain a realistic mass fermion matrix. We have found that the first model 
yields experimentally consistent predictions, where we neglect the Yukawa coulings of the 
first two fermion generations. For the second model we obtain couplings which are so large 
that the model either cannot be treated within the framework of perturbation theory or 
does not give a testable prediction on the top-bottom hierarchy. 

Our approach to predict the top and bottom quark masses, at first sight, looks similar 
to the infrared-fixed-point approach of ref . ||2^ . In the final section of this paper, we will 
discuss this approach within the framework of our concrete 50(10) model and conclude 
that the infrared-fixed-point approach does not always provide us with precise predictions 
on low energy parameters. 
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2 The models 



We denote the hermitean S'O(10)-gamma matrices by Fq , a = 1,---,10. The charge 
conjugation matrix C satisfies C = C^^ , C^^T'^C = — Ta, and the Fn is defined 
as Fii = (— "i)^ n^'LiTa with (Fn)^ = 1. The chiral projection operators are given by 
P± = i(l±Fn). 

In SO{10) GUTs 0, [1^, three generations of quarks and leptons are accommodated 
by three chiral supermultiplets in 16 which we denote by 

^^(16) with = , (1) 

where / runs over the three generations and the spinor index is suppressed. To break 
SO{10) down to SU{3)c x SU{2)i^ x U{1)y, we use the following set of chiral superfields: 

^{„^}(54) , A[„^](45) , 0(16) , 0(16) . (2) 

The two SU{2)i^ doublets which are responsible for the spontaneous symmetry breaking 
(SSB) of SU{2\ X U{1)y down to U{1)em are contained in H^{10). 

2.1 Model I 

For model I, we further introduce a singlet (p which after the SSB of S'O(IO) will mix with 
the right-handed neutrinos so that they will become superheavy. 
The superpotential of model I is given by 

W' = Wy + Wsb + Whs + W},m + Wm , (3) 



where 



^ i,j=i 

WsB = ^0F[,^]0A[„^] + ^Tr ^=^ + ^Tr ^2^, (4) 

Whs = ^ H^S{^p}Hp, W^^^ = J2 9inm'^^'$'P , 
^ 1=1 

Wm = —— H'' + m^Lp^ + m^(p(f) + — S + — A'' , 
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and TiafB] = iiTaTjS — T[^Ta)/2. The superpotential is not the most general one, but 
by virtue of the non-renormahzation theorem, this does not contradict the philosophy 
of the coupling unification by the reduction method (a RG invariant fine tuning is a 
solution of the reduction equation) . Wsb is responsible for the SSB of 50(10) down to 
SU{3)c X SU{2)w X f/(l)y, and this can be achieved without breaking supersymmetry, 
while Whs is responsible for the triplet-doublet splitting of H. The right-handed neutrinos 
obtain a superheavy mass through W^m after the SSB (as announced), and the Yukawa 
couphngs for the leptons and quarks are contained in Wy- We assume that there exists 
a choice of soft supersymmetry breaking terms so that all the vacuum expectation values 
necessary for the desired SSB corresponds to the minimum of the potential. 

Given the supermultiplet content and the superpotential W, we can compute the /? 
functions of the model. The gauge coupling of 50(10) is denoted by g, and our normal- 
ization of the (3 functions is as usual, i.e., dgi/dlufi = jSl;^^ /IQir"^ + 0{g^), where fi is the 
renormalization scale. We find: 



PS 

Phs 

PlNM 

P2NM 

oil) 
P3NM 



7g', 

9t ( 1415^1^ + '^\gHS? + \9mM? - yS-^ ) , 

^ol 48, ,2 1, ,2 1, ,2 1, ,2 77 2 

bi\g^\2 + —\gA\ + -\giNM\ + ^\g2NM\ + ^mNM\ ~ ~^9 



9s{^\9s\ 
5 



l2\gA\' + -\gHs\' - QO9' 



, I ,2 28 I ,2 116 1,2 1 I |2 

^a(16|^0| +—\gs\ +— I^aI +^\gHs\ 



28, 



113, 



gHsiSgr + — gs +M9A + 

5 10 



\9hs\ 

|2 



52/) 
38/) 



45 



.45 2 ni 12 17 2 17, I 

9iNm{—\9(I)\ +y\giNM\ + —\g2NM\ + —\g3NM\ ^ 

/ 45 2 17, ,2 ni i2 17, ,2 45 

g2NM[^\g(j)\ + ^IgiNMl +9|5'2AfM| +y|fi'3iVJ\/| ^ 

gsNMi^lgrl"^ + yls'c/.r + ylS'iiVMl^ + y|5'2iVM|^ + QlgsNMl"^ - y fl'^ ) 



-9 



-9 



(5) 



We have assumed that the Yukawa couplings gu except for gr = (733 vanish. They can 
be included as small perturbations 0. Needless to say that the soft susy breaking terms 
do not alter the /5 functions above. 



^We will clarify later what wc mean by small perturbations. 
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2.2 Model II 

For model II, we introduce a pair of 



0[a/37H(126) and e[a/3^^j.](126) (6) 
instead of the singlet providing us with a possibility of incorporating the Georgi- 



Jarlskog mechanism p5|, |12|- They satisfy the duality conditions 



^ 

©[ai-'-as] (©[oi-'-as]) = ~ ( + ) ^^oi ■■■aio ©[oe-'-aio] {^[ati—aio]) ) (7) 

and G (instead of also will mix with the right-handed neutrinos to make them super- 
heavy. 

The superpotential of model II is given by 

W'' = Wy + Wsb + Whs + Ws + Wgj + W'^'^ + Wm , (8) 

where 

= ^ 0[ai-Q5] 'S'{q5/3i} 0[/3i.../35] + ^ 0[ai-a5] 'S'{a5/3i} 0[/3i-/35] 

WgJ = ^Vl/2CT[„,...,3]Vl/2 0j^^_^, , 

77 1 ^ I — 

^NM = ^ H 9INm'^ C'r[c,,...c,5] 6[oi...„5] ) 

and V]^oi-i_---a5] = -gr (Tai ■ ■ ■ + anti-symmetric permutations ). {Wy , , Whs and 

are given in eq. (4), and Wm contains the 9 — 9 mass term instead of the mass 
term.) The (3 functions are found to be 





= 77 


9', 






= 9t 


2T 63 
( Ulgrl^ + —IgHsl'' + 126|^3iVMr - ) , 






= 94 


48 

53 S'.^p + y P + 63 fifiATAf P + QS\g2NM\'^ + Q^lg^NM 


3 5 1 1 
' + ^\9eA\' 




= 9s 


.84 2^ 2^3 ,2^105 ,2^105 2 


-60/) , 




= 9A 


/-^i i9 28, ,9 116. .9 1. .9 35. .9 
( 16|(7^|2 + —\gs\^ + —\9a\^ + ^\9Hs\^ + y 1^70^1' + 


yl^esl' 



77 2- 
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+f |^70sr-52/) , 

I3hI = ^i/5( 8|f7r|' + y|^7sr + 41^7^1' + ^|^H5r + y|^7e5l' + yl^7e5r- 38/) 

PiNM = 9iNM{Y\9't>\'^ + '^^MaiNAil^ + 71\g2NM\^ + 71\g3NM\'^ + M9gj\'^ 
+yl^7eA|' + 25|(70^|2-^/) , 

+f\geA\' + 25\g^s\^-^g^) , 

PmM = fl'37VM(5|fi'Tp + y |5'</.P + ^llfi-lATA/P + 71|5(2iVMp + 134|5(3ArMp 

nd) / 28 1 |2 ., |2 1| i2 25 1 |2 35 1 ,2 

^es = gesi^lgsl +^\9a\ + ^\9hs\ +y|^eA| + y l^esl 



85, ,2 „ 



2 • 



28 1 
/^§s = ^?e5(yl^?s|' + 4|(/^P + -|f7H5|' + 8|f7Gj|' + 16|(7ijvM|' 

-1/^1 i2 -1/^1 i2 i2 435 I 1 2 35 I 1 2 2 ■ 

+16|fi'27VAf| +16|fi'37VA/| +y|5'eyi| +— l5^5l +yl5'05| -705( 
48 

/^ei = fl'eA( — Ifi-A^ + 8|5(^|^ + 4|5(GJp + 8|5(iArA/P + 8|5(2ArAfr + 8|5(3ArMr 



25 

+l5\geA\' + 25\g^s\' + y l^7es|' - 66^7^ ) , 
Pgj = 9gA 130|(7gj|' + 8|(?i^A/|' + 134|(72iVA//|' + 8\gsNM\^ + y I^^oaI' + 25\g^s\^ - y/ 

Observe the occurrence of large coefficients in the /3 functions above. They are responsible 
for the fact that the model II either cannot be treated in perturbation theory or does not 
give a testable prediction on the top quark mass, as we will see. 



3 Gauge- Yukawa Unification 

The principle of reduction of coupling is to impose as many as possible RG invariant 
constraints which are compatible with renormalizability 0. Such constraints in the space 
of couplings can be expressed in the implicit form as $((71, ■ ■ ■ , (^at) = const., which has 
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to satisfy the partial differential equation 

N fj 

p-v<i> = E = 0, (10) 

where /3j is the f3 function of Qi. In general, there exist, at least locally, indepen- 
dent solutions of (10), and they are equivalent to the solutions of the so-called reduction 
equations 0, 

= A, ^ = 1,■■■,N, (11) 

dg 

where g = qq and (3 = (3q. Since maximally independent RG invariant constraints in 
the (A^ + l)-dimensional space of couplings can be imposed by $j, one could in principle 
express all the couplings in terms of a single coupling, the primary coupling g [^]. This 
possibility is without any doubt attractive, but it can be unrealistic. Therefore, one 
often would like to impose fewer RG invariant constraints, leading to the idea of partial 
reduction [0, |[. 

Here we would like to briefly outline the method Q For the case at hand, it is convenient 
to work with the absolute square of gi, and we define the tilde couplings by 

Oii = — , i = 1, - ■ ■ ,N , 

a 

where a = \g\'^ /At: and = \gi\'^ /At:. We assume that their evolution equations take 
the form 

dt a + , 

—± = -foW + E bi% ajak + ■■■ , (12) 

in perturbation theory, and then we derive from (12) 

«^ = (-i + £)«^-ES«.«^ + E(^r^^r(«), (13) 

where y(\a) , r = 2, ■ ■ -, are power series of and can be computed from the r-th loop 
/3 functions. We then solve the algebraic equations 

(-1 + ^)P^-E-^P.P^ = 0, (14) 



^ Detailed discussions on partial reduction are given in ref. for instance. 



which give the fixed points of (13) at a = 0. We assume that the solutions p^'s have the 
form 



Pi 



for i = 1, ■ ■ ■ , iV' ; Pi > for i = AT' + 1, ■ ■ ■ , iV , 



(15) 



and we regard Ui with i < A^' as small perturbations to the undisturbed system which is 
defined by setting cij with i < N' equal to zero. It is possible to verify at the one-loop 
level the existence of the unique power series solutions 



of the reduction equations (13) to all orders in the undisturbed system (as we will demon- 
strate it in our 5*0(10) model below). These are RG invariant relations among couplings 
that keep formally perturbative renormalizability of the undisturbed system. So in the 
undisturbed system there is only one independent coupling a. 

We emphasize that the more vanishing pj's a solution contains, the less is its predictive 
power in general. We therefore search for predictive solutions in a systematic fashion. 

3.1 Unperturbed system 
(a) Model I 

We find that for model I there exist two independent solutions, A and 5, that have the 
most predictive power, where we have chosen the S'O(IO) gauge coupling as the primary 
coupling: 



P. + Epr(-) 

7^2 ^ 



r-1 



i = N' + 



N 



(16) 





191/204 ~ 0.936 
191/303 ^0.630 



(17) 
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Clearly, the solution B has less predictive power because px = 0. So, we consider below 
only the solution A, in which the coupling a^NM should be regarded as a small perturba- 
tion because Pznm — 0. 

Given this solution, we would like to show next (as promised) that the expansion 
coefficients , i = T,---, '^NM can be uniquely computed in any finite order in pertur- 
bation theory. To this end, wc assume that p\ with r < n — 2 are known, then insert the 
power series ansatz (16) for p\"'^ , i ^ 3NM into the reduction equation (13) and collect 
terms of 0(q;"~^). One finds easily that 



E M,,{n)p^ 



(n) 



known quantities by assumption , i ^ 3NM , (18) 



where 



/ 1141/30 -7n 



M(n) 
















56/81 







283603/9180 - 771 21404/3825 
30467/612 - 7n 152335/4284 

219611/27540 909817/27540 - 7n 
196/405 





31373/1377 


2865/136 
2865/136 

1467/100 


152335/34272 
31373/44064 











28/81 



\ 



(19) 



5351/18360 5351/18360 





-791/810 -7n 

573/68 - 7n 191/24 

191/24 573/68 - 7n J 

So, if det M(n) ^ 0, the coefficients ^ , i = T,---, 2NM can be uniquely calculated. 
We in fact find 



det M{n) = 



110920238635003554634381 3608874567318092545318601 

- + n- 



8522204882112000 



,7571105122486669715209741 
8522204882112000 



— n 



25566614646336000 
, 391617250274453557751579 

284073496070400 
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4 598654192729460650727 g 107001680791 190563 

819127728000 ^ 606761280 

f. 108620968687 ^ 
71^823543 

5508 

7^ for integer n . (20) 
Therefore, there exists a unique power series solution of the form (16) for the solution lA. 



(b) Model II 

According to the principle of reduction of couplings, we search for most predictive 
solutions of (14). Of these solutions, we consider only non-degenerate ones with Pt ^ 0, 
because they are more predictive. We find that there exist three solutions, IIA, IIB and 
lie, and they contain three vanishing p's: 



Pt = 



Ps 



Pes 



PlNM — 



PSNM 



674137 
117840 

108764 
14225 

443 
60 





16361 
4032 



47260924 
4480875 



674137 
2969568 

53321 
159320 

443 
1512 







5.7 
7.7 
7.4 



Pgj 



, Pa 



~ 4.1 



674137 
2969568 

390223 
1433880 

443 
1512 

2356151 , 
5090688 " 





0.2 
0.3 
0.3 



0.5 



, PHS 



4739 
1296 



~ 10.5 

~ 0.2 
~ 0.3 : 
~ 0.3 



Pes 



P2NM — ' 



~ 3. 

673145 
371196 

1584242 
1493625 

16937 
9450 

674137 
2969568 

53321 
159320 

443 
1512 



~ 1.8 
~ 1.1 
~ 1.8 

~ 0.2 
~ 0.3 
~ 0.3 



PeA 



674137 


~ 0.6 


1060560 


300557 


~ 0.6 


512100 


443 
540 


~ 0.8 


1673861 


~ 5.3 


318168 


13811 
51210 


~ 0.3 


77 
81 




~ 1.0 




205226 

36825 — 


= < 


21623024 ^ 
4480875 — 










for < 



IIA 
IIB 

lie 



(21) 



Observe that certain p's for solutions IIB and IIC are so large that the model cannot be 
treated in perturbation theory. The — 5.7 for solution IIA could be within the regime 
of perturbation theory, but as we will see in the next section, that value is so large that 



12 



the predicted value of Mt cannot be distinguished from its infrared value. Therefore, this 
model does not yield a testable prediction on the top quark mass. 

We presented in this section the negative result, too, in some detail to emphasize 
that the existence of a consistent supersymmetric Gauge- Yukawa unified model based on 
50(10) is a nontrivial matter, as we have announced in the introduction. 

3.2 Small perturbations 

The small perturbations caused by nonvanishing cij with i < N', defined in eq. (15) and 
a-sNM in the case of solution lA, enter in such a way that the reduced couplings, i.e., cij 
with i > N', become functions not only of a but also of ttj with i < N'. It turned out 
P that, to investigate such partially reduced systems, it is most convenient to work with 
the partial differential equations, which for solution lA are 

{ PTr + P3~NM^ }ai{a,a) = Pi{a,a) , i 3NM , (22) 

oa oa^NM 

where 

~R - k_L~ Ft - d 

a a a 

These partial differential equations are equivalent to the reduction equations (13), and we 
look for solutions of the form 

&i = P^ + T. fii^3NM) , ^ ^ 3NM , (23) 

r=2 ^ 

where f'f\a^NM) are supposed to be power series of a^NM- This particular type of 
solution can be motivated by requiring that in the limit of vanishing perturbations we 
obtain the undisturbed solutions (16) i.e., fi^\o) = p^p for r > 2. Again it is possible 
to obtain the sufficient conditions for the uniqueness of /j in terms of the lowest order 
coefficients. The proof is similar to that for p^-'^^'s. 

We have computed these corrections up to and including terms of 0{alj^j^,j): 

ar = (163/60 -0.108---a37VM + 0.482---a2^A^ + ---) + --- , 
= (5351/9180 + 0.316---a3iVM + 0.857---a^^A,^ + ■■■) + ■■■ , 
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Figure 1: as versus a^NM-, where the dashed hne is obtained from the analytic expression 



as = (152335/51408 + 0.573 •••a3jvM + 5.7504 •••ag^ji^H )H , 

OiA = (31373/22032 -0.591- --aaiVM -4.832- + - > (24) 

oiHS = (7/81 - 0.00017 •-- + 0.056 -•- a^^A^ + •••) + •• • , 

OiiNM = Oi2NM = ( 191/204- 4.473 •--Q;3ArM + 2.831 ---QigjvA^H )H , 

where - • • indicates higher order terms which can be uniquely computed. In the partially 
reduced theory defined above, we have two independent couplings, a and a^jsiM (along 
with the Yukawa couphngs au , I, J T). 

At the one-loop level eq. (24) defines a line parametrized by asNM in the 7 dimensional 
space of couplings. A numerical analysis shows that this line blows up in the direction 
of as at a finite value of aa^vM- Fig- 1 shows 0:5 as a function of a^NM (the dashed 
line is obtained from the analytic expression (24)). So if we require as to remain within 
the perturbative regime (i.e., gs ^ 2, which means as ^ 8 because q;gut ~ 0.04), the 
Oi3NM should be restricted to be below ~ 0.067. As a consequence, the value of OiT is also 
bounded. To see this, we plot ar as a function of ol^nm in fig- 2, from which we conclude 



(24). 
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Figure 2: ot versus osnm, where the dashed hne is obtained from the analytic expression 



This defines GYU boundary conditions holding at the unification scale Mqut in addition 
to the group theoretic one, ax = c^t = «6 = «r- The value of ay is practically fixed so 
that in the following discussions we may assume that cxt = 163/60 ~ 2.72, which is the 
unperturbed value. 

4 Predictions 

As pointed out, the GYU conditions (25) we have obtained above remain unaffected by 
soft supersymmetry breaking terms, because the (3 functions are not altered by these 
terms. To predict observable parameters from GYU, we apply the renormalization group 
technique [g7[ 0. 

Just below Mqut we would like to obtain the MSSM while requiring that all the 



(24). 



that 



2.714 ^ ar^ 2.736 . 



(25) 
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superpartners are decoupled below the supersymmetry breaking scale Msusy- To simplify 
our numerical analysis we assume a unique threshold Msusy for all the superpartners. 
Then the SM should be spontaneously broken down to SU{3)c x U{1)em due to the 
vacuum expectation value of the scalar component of Ha- We also assume that the low 
energy theory which satisfies the requirement above can be obtained by arranging soft 
supersymmetry breaking terms and the mass parameters in the superpotential (3) in an 
appropriate fashion. 

We shall examine numerically the evolution of the gauge and Yukawa couplings in- 
cluding the two-loop effects, according to their renormalization group equations Q The 
translation of the value of a Yukawa coupling into the corresponding mass value follows 
according to rrii = gi{n) v{jj,)/\/2 , i = t,b,T, where ■mi(fiys are the running masses and 
we use v{Mz) = 246.22 GeV. The pole mass Mi can be calculated from the running one, 
and for the top mass, we use ||22| 



M, = ™,(M,) 1 1 + ii^l^^ + 10.95 i^^r 1 , (26) 

o vr TT 

where we compute v{Mt) from v{Mz)- As for the tau and bottom masses, we assume that 
mT-{fi) and mi,{fi) for /i < Mz satisfy the evolution equation governed by the SU{3)c x 
U{1)em theory with five flavors and use 



on TT 
= m.(M.)[l + ^^MM(i^] , (27) 

TT 

where the couplings with five flavors a3(5f) and Q;EM(5f) are related to and oem by 

«3(5f)(^^) = «3^(^^)-^^' (28) 

«EM(5f)(^^) = «em(Mz)-|^^. (29) 

The corrections in eq. (27) are the SM ones, and in general one should add the MSSM 
corrections too. They could be even large, especially for Mj, ~ ±0(20 — 30%), in the 
case of universal soft supersymmetry breaking terms, while they can be kept small if 



^We take into account the threshold effects by the step function approximation of the (3 functions. 
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these terms are not universal As we will see below, our prediction for mi){Mb) 

without the MSSM corrections fit to the experimental value so that the model favors the 
non-universal soft supersymmetry breaking terms. 
Regarding now 

Mr = 1.777 GeV , Mz = 91.188 GeV , 
a-(M,) = 127.9 + A log ^, 
sin^ OwiMz) = 0.2319 - 3.03 x lO^^T - 8.4 x lO^^T^ , (30) 
T = Mt/[GeV] - 165 , 



as given p3|, we find 



mJMz) = 1.746 GeV , aJMz) = ^ = 8.005 x 10"^ , 

An 

which, together with oem and sin^ 6'w given in (30), we use as the input for the RG 
evolution. In fig. 3, 4 and 5, we show the predictions of model lA on Mf, mfe(Mf,) and 
a^lMz), respectively. Note that the mass values are before the MSSM corrections are 
taken into account. Since mf,{Mh) agrees with the experimental value (4.1 — 4.5) GeV 
23| as we can see from fig. 4, these corrections should be rather small, implying that our 



model favors the non- universal soft supersymmetry breaking terms pl[| . 

Fig. 6 shows that Mf is relatively insensitive to a change in ax- The reason is that 
the predicted values for Mf are not very much far from its infrared value p4|, which we 



define as the value for ctT = Q and is shown in fig. 7. By comparing figs. 6 and 7, we 
see that the predicted values for Mt lie a few GeV below the infrared values. So, if the 
experimental uncertainty can be reduced to less than that order, there will be a chance 
to test the Gauge- Yukawa unified model we have proposed here. Note that the present 
experimental value of Mt is (180 ± 12) [GeV] 

Finally, we would like to comment on the difference of the prediction on Mt with 
and without GYU. As pointed out, the Mf prediction without GYU follows from the 
requirement of the consistent bottom-tau hierarchy. Under the same assumption made 
for the RG analysis above, we wish to find the allowed values of ax- We find that the 
requirement of 4.2 GeV ^ mb{Mi,) ^ 4.5 GeV with MgusY = 500 GeV fixed, for instance, 
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Figure 3: Mt prediction versus Msusy for = 2.717. 




Figure 4: m},{Mi) prediction versus Msusy for olt — 2.717. 
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Figure 5: a^^Mz) prediction versus Msusy for ay = 2.717. 

implies 1.6 ^ ay ~ 7.0, and consequently, 182.9 GeV ^ Mt ^ 191.2 GeV. This should be 
compared with the GYU prediction, 

mb{Mb) = 4.38 GeV , Mt = 187.1 GeV , (31) 

which is fixed under the same assumption. 

5 Comparison with the infrared-fixed-point approach 
and discussions 

The infrared-fixed-point approach is based on the assumption that infrared fixed 
points found in first order in perturbation theory persist in higher orders and that the 
ratio of the compactification scale Ac (or the Planck scale Mp) to Mqut is large enough 
for the ratio of the top Yukawa coupling to the gauge coupling to come very close to its 
infrared value when running from Ac down to Mgut- Since this approach looks similar to 
®In the case of the SM, the infrared fixed point approach loses its meaning at the two loop-level pO|. 
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Figure 6: Mt versus oir with Msusy = 300 GeV (dod-dashed line), 500 GeV (solid line) 
and 1 TeV (dashed hne). 
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Figure 7: The infrared value of Mt versus Msusy (ot = 6). 

ours at first sight, we would like to examine within the framework of our SO (10) model 
whether and how much this picture of infrared-fixed-point behavior is realized. 
To this end, let us first re-capitulate the argument of ref. |^ and recall that 

^(Mgut) 



a (Ac 



(32) 



1 - (7/27r)«(MGUT) ln(Ac/MGUT) 
for one-loop order. For ^(Mgut) = 0.04 and ln(Ac/MGUT) = 5, we obtain a(Ac) — 0.051 . 
Assuming now that a^'s with i ^ T are negligibly small compared with a-p, we derive from 
the reduction equation (13) 



a— = aT{2aT - —) 
da 2 



(33) 



with the solution ||26|, |^ 

""^"^"""^ = 4/11 + (.--4/11)KMgut)MAg)]-/^' ^''^ 

where kt is the value of ckt at Ac- The point in the infrared- fixed-point approach is 

that, since [a(MGUT)/«(Ac)]"^^^^ is small (~ 0.25), the "low-energy" value, q;t(^gut), is 
^ Our normalization of the generators of 5*0(10) differs by a factor of \/2. Our choice corresponds to 
that of the usual SU{5) GUTs. 
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Figure 8: aT(MGUT) (solid line) and ^//^(Mgut) (dashed line) versus a{fi) for kt = 2 
and khs = 2.5. 



insensitive against kt and must be close to its infrared value 11/4 = 2.75. One indeed 
finds that one may vary kt from 1.9 to 4.3 to keep [aT(MGUT)/(ll/4) — 1] ^ 0.1. It 
however should be emphasized that there is no principle to fix in the infrared-fixed- 
point approach. For = 1, which could be realized with the same probability as for 
Kt = 2, we obtain axiMcvT) — 1.91, which is only 70% of the infrared value. 

The more serious problem is the negligibility of other couplings compared to ax- Since 
there exists no reason why the neglected couplings have to be small, they could be large 
and hence comparable to ax, thereby changing the infrared structure very much. In fig. 
8, we plot q;t(^gut) as a function of Q;(/i) with 0.04 < a{fi) < 0.051, where we have 
chosen kt = 2 and khs = 2.5 while neglecting the other couplings in the evolution. As 
we can see from fig. 8, the ar does not approach to 2.75 as a goes to acuT = 0.04 from 
0.051, rather to another fixed point 1.67. (If khs would be zero, the Q!r(MGUT) would 
become 2.51.) 

The observation of ref. that, in spite of the small difference between Mgut and 
Ac, the Yukawa couplings tend to converge to their fixed points very fast thanks to large 
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anomalous dimensions of the matter superfields, is generally correct in one-loop order. 
However, the infrared-fixed-point approach may not always have predictive power, as we 
have explicitly seen above in our concrete model. This is not a specific situation of the 
present SO{10) model, because the factor [a;(MGUT)/tt(Ac)]^^^^ — 0.25 is small enough 



according to the discussion of ref. ||26[. If one insists from the beginning to choose the one- 
loop infrared fixed point which is most predictive, the lowest order prediction is exactly 
the same as that of the reduction of couplings. The difference appears in the next order, 
because the reduction solution, except for the lowest order, is not a fixed point solution 
in general. 

If the low-energy prediction for either approach is viable as experienced in some 
cases including the one discussed in this paper, it might indicate some unknown non- 
perturbative mechanism of unification of couplings such as the dynamical unification of 
couplings which we suggested in introduction. In either case, the non-perturbative inves- 
tigation on asymptotically-nonfree, non-abelian gauge theories will be an important issue 
in future works. 
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